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Abstract 

Understanding the power spectrum of the magnetization noise is a long standing 
problem. While earlier work considered superposition of 'elementary' jumps, without 
reference to the underlying physics, recent approaches relate the properties of the 
noise with the critical dynamics of domain walls. In particular, a new derivation of 
the power spectrum exponent has been proposed for the random-field Ising model. 
We apply this approach to experimental data, showing its validity and limitations. 
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After almost a century since its discovery in 1909, the magnetization noise 
produced by the intermittent motion of domain walls, i.e. the Barkhausen 
noise, still represents an intriguing scientific challenge from the theoretical 
point of view. Considering the long-time and vast production of experimental 
and theoretical papers of the past, it is quite surprising that only recently 
an exhaustive comprehension of the noise properties has been achieved. The 
introduction of methods of statistical mechanics, in fact, made possible a re- 
liable description of the intrinsic complexity of magnetization processes. In 
particular, the power law exhibited by the Barkhausen signal amplitude to- 
gether with the avalanche size and duration has been explained in terms of 
an underling critical point. Its true nature is still under debate, as two main 
different approaches have been proposed: the zero temperature random-field 
Ising model (RFIM) [1] , where criticality is set by the amount of disorder, and 
interface model, where a domain wall moves through a disordered medium 
and criticality is due to the depinning transition of the wall [2-4]. Exploit- 
ing the effects of long range dipolar magnetostatic fields and of domain wall 
elastic tension on the depinning transition, we have recently shown that two 
distinct universality classes exist, with different critical exponents and cutoff 
scaling dependence on the demagnetizing factor due to sample geometry [4]. 
These results have been confirmed on two polycrystalline and amorphous sets 
of materials, supporting a domain wall theory for the Barkhausen effect. 

Despite these significant results, a proper description of the shape of the 
power spectrum noise is still not available. Significantly, about 80% of the 
Barkhausen literature have been devoted to this problem. Earlier approaches 
considering a description of the power spectrum shape as a superposition of 
elementary independent events (see for instance [5,6]) appeared to be unre- 
lated to any microscopic mechanism, thus not clarifying the true origin of the 
magnetization process. Even attempts to link the power spectral exponent to 
the critical exponent of size distribution by a simple scaling relation [7-9] ap- 
pears to be quite unsatisfactory and not confirmed in general by experiments. 
Experimental noise in fact show a quite complex pattern: the high frequency 
part often follows a power law with exponents between 1.5 and 2, even if 
cases have been reported where the power law extension is very limited and a 
more complex pattern results [10]. The low frequency part displays a marked 
peak, those position strongly depends on the magnetization rate, and a 
dependence at lower frequencies, with 0.5 < a < 1. Such a complicated pat- 
tern does not actually reveal all the complexity of the underlying dynamics of 
Barkhausen avalanches. Considering high order moments of the signal, differ- 
ent frequency bands appear to be very strong coupled, so that also high order 
power spectra display power law dependences [11,12]. In addition, time asym- 
metries of third-order voltage correlations are found in amorphous samples, 
showing as high frequency events precede on average the low frequency ones 
[11]. All these properties are not currently explained by any of the existing 
models. 
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A derivation of the power spectrum exponent from a scaling analysis, at least 
for some simple models, can help for a better comprehension of the avalanche 
dynamics. A decisive step in this direction has been performed recently by 
Kuntz and Sethna [13], who derived the power spectrum exponent in the zero 
temperature RFIM, and corrected the earlier theoretical estimations [7-9]. 
Most part of their derivation can be applied to different models, and also to 
the analysis of real data. In this paper, we apply some of their results to various 
sets of experimental data, confirming the general validity of this approach for 
simple cases, and showing where it must be improved for a better description 
of the results. 

Let us consider some of the relations of Ref. [13] that we will use in our 
analysis. The key scaling relation connects the average avalanche size (s(T)) 
with its duration T, that is {s{T)) ~ j^'^/'^^z ^^ere the exponents cr, v and z 
are defined in Ref. [8] . When the exponent of the avalanche size distribution 
r is less than 2, as usually in experiments [4], the high frequency tail of power 
spectrum is calculated to scale with exponent 1/avz. This central result is 
based on the existence of a couple of scaling relations regarding the avalanche 
shape. The first one states that the average avalanche shape should scale in a 
universal way, so that 

y(T,t) = r^/-^-V,,„pe(t/T) (1) 

where V is the signal voltage, t is the time and fshape{t/T) is a universal 
scaling function, having the approximated shape of an inverted parabola for 
the RFIM. The second relation analyzes the fluctuations of avalanche sizes 
considering the probabihty P{V\s) of the occurrence of voltage V inside an 
avalanche of size s. This probability scales as: 

P(V\S) = V-'Ultage{Vs-''^-') (2) 

where fvoUage is another universal scaling function. With relations 1 and 2, 
the power spectrum exponent is obtained calculating the time-time correla- 
tion function in the case of adiabatically increase of the applied field and of a 
complete separation of avalanches in time, thus avoiding any avalanche corre- 
lation. All these considerations will be helpful to understand the experimental 
results. 

We consider two kind of samples, belonging to different universality classes as 
pointed out in Ref. [4]. An as-cast FeQ4Co2iBi^ amorphous alloy, measured 

under moderate tensile stress {a ~ 20MPa), and an Fe-Si 7.8 wt.% strip 
(30 cm X 0.5 cm x 60 /im) produced by plan flow casting, having grains of 
average dimension of 25 ixm. The amorphous ribbon follows the universality 
class where the surface tension of the wall dominates the domain dynamics 
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(short range class), where Xjavz ~ 1.77 [4]. Fig. 1 shows the comparison of 
the power spectrum with the average size distribution (s(l/T)) as a func- 
tion of the inverse of avalanche durations. The agreement with the theoretical 
prediction is fairly good over an extended time range: at high avalanche du- 
rations (small frequencies), the time correlation between avalanches becomes 
relevant and the theoretical analysis is no longer valid. The inset of Fig. 1 show 
the same comparison in the case of the FcSi sample. This material falls in the 
universality class where long range magnetostatic fields dominates the domain 
dynamics, giving l/avz ~ 2 [3]. Also in this case, the agreement is fairly good, 
but for a smaller high frequency range. The precise reason for this fact is not 
clear, even though a visual inspection of the time signals of both materials can 
justify this result: in the amorphous alloy, the avalanches are well separated 
in time (see [14]), while in the FeSi alloy the separation is much less defined 
(see [15]). This is also confirmed by the fact that only in the latter material the 
critical exponents of size and duration distributions strongly depend on the 
applied field rate. This means that avalanche correlations, and thus time-time 
correlations, are significantly different, giving different spectral contributions. 
The results for the amorphous alloy are confirmed by the scaling of the av- 
erage avalanche shape V{T,t) (fig. 2) and of the probability P{y\s) (fig. 3). 
In both figures the theoretical value 1/avz ~ 1.77 is used. Interestingly, the 
universal scaling function fshape of eq. 1 is not an inverted parabola as for the 
RFIM, but it shows a marked temporal asymmetry. This is compatible with 
the results of Ref. [11] concerning high order spectra: the high frequency signal 
components precede the low frequency ones, so that an avalanche of a given 
size and duration starts with a fast signal ramp and relaxes at longer times. 
Interestingly enough, despite this time asymmetry, the predictions of Ref. [13] 
are confirmed, suggesting that the scaling properties are more important than 
the exact shapes of average avalanches. This conclusion strongly contradicts 
with the basic assumption often reported in the literature [5,6], where a distri- 
bution of 'elementary' avalanches with a pre-defined shape (often exponential) 
is summed up to calculate the power spectrum. 

We must add that the avalanches of the FcSi alloy does not show such a 
nice scaling. In particular, short and long avalanches have markedly different 
shapes, as the former is approximately an inverted parabola, while the latter 
show a flat central region. Also the scaling of P{V\s) is not perfectly compat- 
ible with the theoretical exponent l/auz = 2. All these features could explain 
the limited agreement between the power spectrum and (s(T)), and surely 
need a more extensive analysis, taking into account avalanche correlations 
and the dominant role of demagnetizing fields. 

From the analysis of experimental results shown above, one may argue that 
materials belonging to the short range class also exhibit power spectra scaling 
as l/auz. As a matter of fact, the application of larger applied tensile stresses 
on the amorphous material does not change the universality class [14], but 
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reveals a more complex behavior. In particular, S{l/u) and (s(T)) do not 
longer scale in a similar way at high stresses, even if V{T,t) and P{V, s) still 
rescale approximately with 1/avz = 1.77. This behavior reflects the change 
of the avalanche shape as shown in Ref. [14]: the cutoff of avalanche duration 
decreases while keeping the size distribution invariant. As the scaling range of 
duration distribution gets shorter, we expect that the longest avalanches, close 
to the cutoff value, are increasingly more effective in the time-time correlation 
of the signal, and thus to the power spectra. Their frequency content is not 
obviously taken into account in the scaling calculations of [13]. 

As clearly shown above, many different aspects can enter in the definition 
of the scaling properties of the Barkhausen signal. As already pointed out 
[11,12], a more complex analysis is required to evaluate in detail the statistics 
of avalanches. In this respect, the analysis of Ref. [13] not only gives a new 
approach to the long standing problem of power spectra shape, but introduces 
precise statistical tools helpful both to the experimental analysis and to the 
theoretical description of magnetization processes. Within some limitations to 
be further investigated, the results described above are fully compatible with 
our interface models [2-4] . This is quite surprising, considering that this model 
considers the motion of a single domain wall with strong simplifications of the 
fields relevant to the dynamics (especially of magnetostatic fields). In fact, a 
real material typically shows a complex pattern of multiple domains with no 
easily predictable field configuration. Here temporal and spatial correlations 
are important, resulting in a very complex behavior of time-time correlation 
and thus of power spectra. 
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Fig. 1. Comparison of the power spectrum ^(a;) with the average avalanche size 
s(l/T) as a function of inverse of avalanche duration T for an FeQiCo2iBi^ amor- 
phous ribbon and an FeSi 7.8 wt.% polycrystalline alloy (inset). The theoretical 
prediction of Ref. [13] is also shown, with 1/auz equal to 1.77 and 2, respectively, 
as given by interface model [4]. 




Fig. 2. Average avalanche shape (eq. 1) for the Fe64Co2i-Bi5 amorphous ribbon 
using the scaling exponent 1/avz equal to 1.77. The full line is the average of the 
curves and the dotted line is a symmetric parabola. 
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Fig. 3. Distribution of avalanche voltages at fixed voltage (eq. 2) for the Fe^Co^iBi^ 
amorphous ribbon using the scaling exponent 1/avz equal to 1.77. 
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